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Abstract 

We investigate the mass, production and branching ratios of a 125 GeV Higgs 
in models with Dirac gaugino masses. We give a discussion of naturalness, and de- 
scribe how deviations from the Standard Model in the key Higgs search channels can 
be simply obtained. We then perform parameter scans using a SARAH package up- 
grade, which produces SPheno code that calculates all relevant quantities, including 
electroweak precision and flavour constraint data, to a level of accuracy previously 
impossible for this class of models. We study three different variations on the minimal 
Dirac gaugino extension of the (N)MSSM. 



1 Introduction 



The LHC experiments have claimed the discovery of a new particle at 125 ± 1 GeV [l| 
[2]. Its production and decays make it a good candidate for the Higgs boson. However, 
more experimental data is needed for a precise determination of its quantum numbers 
and interactions to know if the Higgs sector is the Standard Model one, or an extended 
version, for instance as required in supersymmetric models. A clear indication of 
non-minimality would be a significant excess or deficit in at least one decay channel, 
but the Higgs mass itself can also be regarded as such, since it is somewhat high 
for the most natural version of the MSSM. This is one motivation for this work, 
in which we study the main properties of such a Higgs boson in models with Dirac 
gaugino masses [3 -31 



These have numerous virtues compared to their Majorana 



counterparts, not least that they allow for increased naturalness [6,21 



36 ; but also 



is particularly important in the light of recent LHC SUSY searches [35 
that the direct production of gluinos is suppressed, loosening the bounds from direct 



32-34 which 
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LHC searches 27,29 ; they allow relaxation of flavour constraints such as due to 
Br(i? — )• 57) [37] ; they preserve R-symmetry (allowing for simpler supersymmetry- 
breaking models) and can be motivated from higher-dimensional theories as being 
derived from an = 2 supersymmetry in the gauge sector at some scale. 

To give gauginos Dirac masses, we must add new chiral superfields in the adjoint 
representation of each gauge group: a sing let S = S" + V2exs + • • • , an SU{2) triplet 
T = Ea=i,2,3 T^"^ = + V2exT^ + • • • and an SU{3) octet Og = Ea=i..8 O^"^ + 
V^ex'o+- • • . Of these, the triplet and singlet may have new renormalisable couplings 
with the Higgs which allow more possibilities to obtain the desired mass range than 
in the MSSM. More precisely, the lightest Higgs mass in these models is determined 
by opposite competing effects. On one hand, the presence of new couplings in the 
extended scalar sector leads to an enhancement of this mass by allowing new contribu- 
tions to the quartic Higgs coupling. On the other hand, the supersoft operators that 
include the Dirac mass induce new D-term couplings, which increase Higgs mixing 
and thus tend to reduce the lightest Higgs mass. However, this is only potentially 
problematic if the triplet scalar soft mass is small; but we shall demonstrate in section 
[4]- along with a general discussion of naturalness in Dirac gaugino models - that it 
may be naturally large enough to avoid this problem. Higgs mixing involving the sin- 
glet induced by a Dirac Bino mass, however, then presents an intriguing opportunity: 
it allows the decays of the Higgs to b quarks and r leptons to be suppressed while 
preserving the rate to Ws and Zs, and an enhancement of the diphoton rate. We give 
an explanation of this and a discussion of the Higgs production and decay rates in 
section [H 

Previous attempts to quantitatively study the Higgs sectors of Dirac gaugino mod- 
els have been hampered by the lack of numerical tools to do so; until now only one- 
loop effective potential calculations of the Higgs mass have been possible. This is in 
contrast to the MSSM, where the leading corrections are known to three-loop order. 
However, with an upgrade of the SARAH package described in section |3j it is now possi- 
ble to study models with arbitrary gauge groups and matter content: it can generate 
SPheno code which calculates all one-loop pole masses and tadpoles, which allows 
a much more accurate determination of the Higgs mass. We implement a minimal 
Dirac-gaugino extension of the (N)MSSM which lends itself to four particularly inter- 
esting sub-classes of models, which we review in section [2] One such class of models 
is the "MSSM in disguise," where all the new scalars are too heavy to observe or mix 
substantially with the Higgs; this is a good toy scenario to use to test of the code, 
and we do just that. 

The SPheno code produced by SARAH can also calculate the branching ratios and 
production cross-sections of the Higgs, as well as electroweak precision observables 
such as A/9 and flavour constraints such as Br(i? — t- 57). We take full advantage of 
the latter in investigating the "MSSM without fi term" ^ in section [6| comparing the 
results for Ap with approximations given in |7]. Unfortunately those constraints in 
addition to those on chargino masses and the Higgs mass yield that model problematic. 
However, we propose a new model, which we call "dynamica n models," in which the 
singlet obtains a substantial expectation value - we show in model scans in section [7| 
that this not only alleviates all of the problems of the MSSM without fj, term, but also 
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allows for Higgs mixing and thus interesting deviations in the Higgs production rates 
and branching ratios that may better fit the current data than the standard model. 

2 One model into four 

Adding Dirac gaugino masses to the MSSM introduces many extra parameters: not 
just the Dirac masses themselves, but also new superpotential couplings and soft 
terms. By making certain assumptions about the new parameters we can then arrive 
at different limits of the model with different phenomenology; we shall consider four 
such limits. 

The most general renormalisablc superpotential that we can write down is 

W = YuUqH^-YjqHa-Y.elHa + fiHu-tld 
+A5SHu-Hd + 2ATHd-THu 

+LS + + ^S^ + MTtr(TT) + Motr(OO) 

+W2 (2.1) 

where 

W2 = A5TStr(TT) + AsoStr(OO) + ^tr(OOO). (2.2) 

The usual scalar soft terms are 

-A/:|^fsV°^* =[TuUqHu - TddqHd - TeelHa + h.c] 

+ mljH^l'' + mljHdl'' + • Ha + h.c] 

+ q\ml)% + e{ml){ui + (m^)^dy + ^(mf )^Z^ + &{rr^e^ (2.3) 
and there are soft terms involving the adjoint scalars 
-A/:sfo1ntf* = {tsS + h.c) 

+m||5p + \bs{S'^ + h.c) + 2m|.tr(r1'T) + (Brtr(Tr) + h.c) 

+ [As\sSH^ ■ Hd + 2AT\THd ■ THu + ^/^^^^^ + h.c] 
+2m^tr(Oto) + {Boti{00) + h.c) (2.4) 
with the definition Hu ■ = H+H~ - H^H^. Similarly there are the A-terms for 

One limit of the above model would be to allow all parameters, including Ma- 

jorana gaugino masses, to be significant and non- vanishing. Such models may have 
virtues due to the extra Higgs couplings and extra gaugino states (the charginos could 
contribute, for example, to enhancing the Higgs-to-diphoton decay channel) but we 
shall leave the exploration of this to future work. Instead, we shall explore models 
where the gaugino masses are entirely Dirac, taking as motivation the possibility of 
preserving R-symmetry in some sector of the theory (but not in all: it must ultimately 
be broken in any case). Of particular interest to us are: 
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• MSSM in disguise: here we shall allow a /x-term, and assume that the only source 
of R-symmetry violation arises in the supersymmetry-breaking sector, but per- 
mit only a i?^ term. This assumption will be preserved by renormalisation group 
running and so it can be justified by high-energy boundary conditions; non-zero 
A-terms would lead to Majorana gaugino masses. Since it is the "MSSM in 
disguise" we shall take the scalar singlet and triplet to be very massive (several 
TeV) . We perform some scans over models of this type in section [s] 

• MSSM without /i term (^SSM): this is the scenario of |7|, similar to the MSSM 
in disguise but taking = 0. Here we shall insist that the singlet vev is small, 
so the chargino mass must be generated by the coupling At- We investigate this 
scenario in section |6l 

• Dynamical /j, models: in this scenario, we again take /i = but allow non-zero 
Bs, leading to a substantial non-zero expectation value for the singlet. In this 
way, the vev and the coupling A5 lead to an effective //-term. Models of this 
type are very natural and interesting from the point of view of Higgs mixing: 
we perform scans over these models in section [7j 

• Dynamical /x and Bfi models: this is the scenario of [21] where we allow a non- 
zero K, breaking R-symmetry in the visible sector, but allowing /x and to be 
generated via a non-zero singlet vev - so we can set all R-symmetry-breaking 
parameters in the supersymmetry breaking sector to zero. It is somewhat sim- 
ilar to the NMSSM, but the Dirac masses lead to some interesting differences. 
Models of this type can be very natural, but we leave scans of their parameter 
space to future work. 

We review the tree- level properties of the generic case in appendix [Aj see also [2l| . 
However, common to all of the above scenarios is the assumption that Dirac gaugino 
masses dominate over Majorana ones. For this to be natural, A-terms must also be 
small, so in our searches we shall set (unless otherwise stated) W2 = L = Ms = Mt = 
Mo = As = A^ = At = 0. 



Some notation 

We now introduce some notation relevant for the following: we redefine the singlet 
and triplet scalars in terms of real components S = -^{vs + sr + isj),T^ = :^(t'T + 
Tji + iTj) and have an "effective /i-term" fl = fi+ -^{vsXs + VT)yT)- The expectation 
values vs,vt are associated with new non-trivial potential minimisation conditions, 
which we give in equation (A. 7). The scalars sr, Tr mix with the Higgs fields, with a 
4x4 mass matrix given in equation (A.12); the mixing will be important in section 

El 
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3 Numerical Setup 



3.1 Implementation in SARAH and SPheno 



For our numerical studies we have extended the Mathematica package SARAH [38-41 
to support Dirac gauginos and implemented our model Afterwards we used the 
possibility of SARAH to create source code for SPheno |42[|43| for the given model. 
The obtained code was compiled with SPheno 3.1.2 and provides a fully functional 
spectrum calculator: the entire mass spectrum is calculated at one-loop using the 
full dependence on the momentum of the external particle. For a detailed discus- 
sion of the calculation of the loop corrected mass spectrum using SARAH and SPheno 
for extensions of the MSSM we refer the interested reader to Ref. |44] . In addition, 
the produced SPheno version includes routines to calculate the decay widths and 
branching ratios for all SM superpartners and Higgs fields. In general, these are pure 
tree-level calculation. However, in the case of the Higgs particles the loop induced 
decays into two photons and two gluons are included. This calculation is comparable 
to the analytical results given in sec. [5j but also the NLO corrections due to quarks 
and squarks are added as given in Ref. [45]. Similarly, for the Higgs decays into 
quarks the dominant QCD corrections due to the gluon are taken into account (46]. 
Finally, the SPheno version for Dirac gauginos also calculates the observables 6 — )• 57 
and Ap with the same precision as done in the MSSM by SPheno 3.1.1 including all 
contributions due to the new states present in the considered model. The details of 
these calculations can be found in Ref. |43| and references therein. 

Since there are four non-trivial vacuum minimisation conditions (given in equa- 



tions (A. 5), (A. 6) and (A. 7)) we must use these to eliminate four parameters. From 



a top-down perspective, it would be preferable to specify the soft masses ms,mT as 



was done in 21 , and derive from that vs,vt- However, the equations are non-linear 
in these, and so it is preferable to instead take vs,VT,tanf3, fj, as input parameters 
in the code, and treat m'g,m'^,m'\j^,m?^^ as output parameters (calculated includ- 
ing one-loop tadpoles). As further inputs for our studies we use for our studies the 
soft-breaking terms as well as the superpotential couplings at the SUSY scale: the 
standard model gauge and Yukawa couplings are calculated at MsusY using two-loop 



standard model RGEs from Mz 47 



Finally we note that, in the absence of the theoretical calculation, the code cannot 
include the two-loop corrections to the Higgs mass involving (Dirac) gluino masses 
that can be important in the MSSM and NMSSM. In those cases, the Higgs mass 
is usually increased by 2 to 4 GeV. Throughout we shall be conservative and allow 
a variation of ±4 GeV for the mass of the Higgs in the scans, but we expect that 
typically the shift will be positive. 



"'^The support of Dirac gauginos as well as the used model files will become public with the next major 
upgrade to version 3.2.0 of SARAH. 
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3.2 Comparison with effective potential 



It is possible to calculate an approximation for the Higgs mass via the effective po- 
tential technique. In the decoupling limit of large Bfx, this yields 



y2 



+ 



3 



log^ + ^™^|l 



mr 



cot^ /3 



XlCp + X2s'p + 2(A3 + A4 + A5)c|s| + 4.{XqCp + A7S^)S/3C/3 



tan (9— >oo 



3 



, mr mr 

Mog^V^ 

mt 



(3.1) 



where the first line is the tree- level mass, the second line is the usual contribution from 
stops with mi_^ ^ the physical masses and At set to zero; the A, are the coefficients of 
the terms in the most general CP-conserving gauge-invariant potential up to quartic 

\HyY^ operator) about zero vev 48 . This 



order (e.g. A2 is the coefficient of the \Hu\'^ operator) about zero vev 
is a good approximation when the Higgs vev is small compared to the energy scales 
being integrated out; in this case, we shall consider integrating out the adjoint scalars, 
which is appropriate for the MSSM in disguise. We give the full expression for the 
coefficients Aj in appendix [B| but an interesting limit is to consider Bg = Bt = and 
neglect v,vs and the Dirac masses rnDf, then the contribution from the singlet and 



triplet scalars is 21 



+ 



"4 
4X^X1 



32n'X2 D2A| log '-^ + {9^- ^g^X^r + lOA^) log 



nin 



AglX^T + 2(A| 



m. 



■log- 



mt 



mr 



■log- 



mn 



(m| 



m^) 



2Ac^A^ ~t~ SA^-i 



log- 



mt 



(3.2) 



where in the last line we show the limit that the scalar masses become equal. 

By performing a scan over parameters we can compare this approximation with 
the more accurate results from SPheno using the code produced by SARAH; we show 
the results with the choice of tan/3 = 50, mD2 = 600 GeV, first two generation 
sfermion mass squareds of 4 x 10^ (GeV)^, third generation sfermion mass squareds 
4 X 10^ (GeV)2 and scanning over moi e [-600, 600] GeV, fi G [-750, 750] GeV, Bfi G 
[5000, 10^] (GeV)2, At G [0, 1] while adjusting Xs to keep mh = 125 ± 4 GeV in figure 
|3.2[ The expectation values vs,vt were set by the tree level minimisation equation 
with m|,m^ = 2.5 x 10^ (GeV)^; this resulted in values close to this for m^, while 
the resulting one-loop adjusted values for m| varied between 10^ and 10^*^ (GeV)^. 
As can be seen from the figure, there is good agreement between the two, although 
of course the approximate effective potential calculation exhibits a wider variation of 
masses; that there is no apparent correlation is unsurprising, essentially reflecting the 
error margin in the effective potential method. 



6 



122 124 126 128 

Approx m /, (GeV) 

Figure 1: Comparison of effective potential method (labelled "approx m^") with a parame- 
ter scan using SPheno code produced by SARAH. Parameters are given in the text; a searched 
was performed using the SPheno code, restricting to nth = 125 ±4 GeV, and from the same 
input parameters for each model the one-loop Higgs mass is shown (note that this has a 
larger range). Only models with relatively low couplings A^, Aj- (roughly A| + A|- < 1.2). 



4 Dirac gauginos and Natural SUSY 

One of the increasingly attractive features of Dirac gaugino models is that the su- 
persoft operators allow for increased naturalness (6| [2H[32[|36] ; they do not appear in 
the RGEs for the soft masses and so they only affect the stop mass via a UV-finite 
correction, allowing for heavy gluinos with light stops. On the other hand, the Higgs 
potential is corrected at tree-level by two competing effects: one is the enhancement 
of the Higgs mass (at low tan /3) by the new couplings Ag, A-p, clearly evident in equa- 



tion (3.1); the other is a reduction in the effective D-term Higgs quartic coupling due 
to the Dirac mass terms; if we ignore the superpotential couplings and integrate out 
the heavy scalars then |6lpi] 



g^ + g^ If ml + \Ms\^ + Bs 2, ml + \Mt\'^ + Bt 



i I J ■ I ^ I ■ 2 

ml + \Ms\^ + Bs + 4|mmP^^ m^, + |MtP + Bt + 4|mD2P^^ 

(4.1 



This latter effect manifests itself as mixing terms in the Higgs mass matrix (A. 12) 
Ahs,Aht; provided that the other soft masses are large enough the suppression can 
be avoided. We may then ask how large these soft masses can be while preserving 
naturalness, to which we consider the ratiative corrections to m?r : 

" 16^^^^^""^ ^ 6A^m^) log (4.2) 
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where A is the UV cutoff of the theory. Using A = then we have 



ms <TeV 



1 \ f\ogA/TeVY [ ^ 



1/2 / A-l\ 



2 



\sj \ 3 J V20% 
„.,..v(^)(^)--(-)--. ..3, 

The latter value is particularly important, because we also have the constraint that the 
triplet scalar expectation value must be small in order to avoid a large p parameter. 
We have (see appendix) 

^2 

- I 2 ^1,^ 12^ R V2V2^92niD2C2p + \/2/iAT)' < 8 X 10"^ (4.4) 

which is satisfied for any value of mj:,2 for mj^ ^1-4 TeV; but interestingly is also 
satisfied for any value of rriT if mD2 ^ 1-4 TeV. 

An interesting choice for the parameters \s,\t would be for them to take their 
N = 2 supersymmetric values at some scale, Xs = gy/V^jXr = g2/V^, although 
corrections due to the running would break this exact relation. However, even if the 
N = 2 scale is intermediate (such as 10^^ GeV) assuming a desert we would still 
find Xs — 0.2, Xt — 0.4 at the low energy scale. In this case, in order to preserve 
naturalness we would require 

In order to satisfy the tree-level Ap constraint, we would either need to work at small 
C2I3, or take mci2 ^ 1-4 TeV and ensure that the ensuing Higgs mixing allows a large 
enough Higgs mass. We leave exploring this interesting possibility to future work: 
in this paper we shall take large values of my and small values of At to keep the 
loop-level corrections to Ap small. 

One final naturalness-related issue in these models is that the requirement of 
small ^-terms typically diminishes stop mixing; in the case of large tan /? the mixing 



is almost eliminated. Hence the contribution from the stops to Ap is 49 -51 



^ stops ^ 3ae^ 2 2 , 



ap,m mf 4 X 10^4 I^SOOGeVy 



(4.6) 



(where Fo{a, b) = a + b — log |) which is similar to the experimental value. In 
"natural SUSY" MSSM models the stops are lighter than about 600 GeV, and so the 
stops by themselves will not be able to lift the Higgs mass to 125 GeV, but have a large 
impact on the electroweak precision corrections. On the other hand, in the context 
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of ASUSY 51 the model can remain natural for stops as heavy as 1.5 TeV 51,52 



because the relative contribution of stops to the Higgs mass compared to the tree-level 
effect is small. Clearly A5 in our models is the same as the A in ASUSY. Hence the 
simplest natural scenario is to take small tan /3, small At, and As ~ 1 with ~ few 
TeV and ms ^ TeV. The Dirac gaugino masses mD^-,'mD2 can be naturally large, but 
the Dirac Bino mass will be bounded above by the amount of mixing that it induces 
between the singlet and lightest Higgs (thus reducing the Higgs mass) to be a few 
hundred GeV. As a result of this, an amusing feature is that natural Dirac gauginos 
will lead to a Majorana neutralino, since there will be non-negligible mixing between 
the Bino and the neutral Higgsinos. 

It is therefore worthwhile giving an example of such a "natural" model. In the 
context of "dynamical models we take A^ = —l.O,B^ = 8.0 x 10'^ (GeV)^,?;^ = 
170 GeV and mm = -170 GeY,mD2 = mos = 10^ GeV, mo = 3 TeV, ts = -1.5 x 
10^(GeV)3, Bs = -5x 10^ (GeV)^, vt = 0.465 GeV, = 2.2 x 10^ (GeV)^ while the 
first two generations have soft mass squareds of 4 x lO"^ (GeV)^, with third generation 
soft masses at (500 GeV)^, which leads to light neutral Higgs masses (at one loop) of 
121, 362, 429, 518 GeV, a light pseudoscalar Higgs mass of 438, a light charged Higgs 
mass of 412 GeV, neutralinos of masses 82, 134, 217, 267, 1056, 1056 GeV and charginos 
of masses 121 GeV and a TeV. We then find Ap = 8 x 10"^, BR(6 57) = 3.2 x 10""^. 
We shall now turn to a discussion of the Higgs signatures at the LHC of natural Dirac 
gaugino models. 



5 Higgs production and branching ratios 

It is now important to consider the production cross-sections and branching ratios of 
the Higgs. In our Dirac gaugino models there is a singlet scalar which may mix with 
the lightest Higgs state, so we shall consider the branching ratios into each channel 
taking into account the mixing, providing some approximate expressions and then 
comparing them to the output of the SPheno code created by SARAH. 
We shall use the standard definitions 



Ri = 



BR{h ii) 



BRsuih — )■ a 
(t{pp — )■ h) 



Ri 



(5.1) 



where An is the amplitude. 

To take Higgs mixing into account, consider the rotation of the states /i, sr via 
a matrix S so that 



h ^ 




f hi 


H 






SR J 







(5.2) 



We shall assume throughout that the lightest Higgs field hi has mass 125 GeV, i.e. 
there is no additional lighter singlet. In this notation, we can then calculate how 
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the production and decay channels are modified. At 125 GeV the production cross- 
sections (as Hsted on the CERN yellow pages [53]) are 



o'suipp — > h) =19.5;|;J^'^^°pb gluon fusion 

-|-1.578^g;^^°pb vector boson fusion 
+0.6966+^-JJ°pb WH process 
+0.3943+^-5%pb ZH process 
+0.1302+}^;^|pb ttH process 

Of the initial eigenstates, only h couples at tree level to the vector bosons. The 
coupling to gluons and in the ttH process is proportional to the top quark coupling 
squared, so 

r{gg hi) ~ hiti oc + cot l3\^. (5.3) 

Hence we can write 

^{gg hi) (19.5 + 0.1302)|5ii + 521 cot /3|2 + (1.578 + 0.6966 + 0.3943)|5ii| 



' Tsuigg^hi) 22.2991 

=0.88|5ii + 521 cot /3|2 + 0.12|5iip. (5.4) 

We can use the same approach for the decay branching ratios: 

Rg, Rc oc htt ~|5ii + 521 cot 



Rb oc T{hi bb) ~ 
Rr OC T{hi —7- fr) ~ 
Rw oc r(/ii ^ WW*) oc 
Rz oc T{hi ZZ*) oc 



5ii — 521 tan /3|^ 
5ii — 521 tan 
5iip 

5ii|'. (5.5) 



Since the photon couples at one loop to the Higgs, and the singlet couples to charged 
fields, the expression for the coupling to the photon will be more complicated. At 
nih = 125 GeV, the standard model Higgs branching ratio (as listed on the CERN 
yellow pages [53]) is dominated by 

BRsM{hi bb) =5.77 x 10"^ 
BRsM{hi WW*) =2.15 X 10"^ 
BRsMihi ^ gg) =8.6 x 10"^ 
BRsuihi tt) =6.3 X 10"^ 
BRsAiihi cc) =2.91 X 10"^ 
BRsAiihi ZZ*) =2.6 X 10"^ (5.6) 

So we can write 

i?x(l — Br(/i — > invisible)) 
" ' ErRYBrsMih^Y) 

-Rx(l - Br(/i invisible)) 

~ ^0.640|5ii-52itan/3|2 + 0.241|5ii|2 + 0.115|5ii + 52icot/3|2' ^ " ^ 
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To approximate the diphoton channel, we require the expression 



9hiVV ^2 



m 



QpAi{Tv) + 



2gh,ffNcQ) 



V 



ruf 



Al/2{Tf) + 



ghiss^cQs 



Mrs) 



(5.8) 



where the functions Ag where s is the spin of the field in the loop are standard and 
given in 



45 



Tx ^ 4m^/m^^. For fermions and scalars they are well approximated 
by the large mass limits of 4/3 and 1/3 respectively, but for the W boson and top 
quark the values are Ai{tw) — —8.32, (4/3)Ai/2(''"t) — 1.84. If we consider that the 
singlet couples to some charged Dirac fermion we can write for the couplings ghiU 



2M, 



dhiVV 
Qhiti 
9hibb 



V 

mt 

V 

nib 
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(5ii + cot/352i) 
{Sii - tan/3S'2i) 



(5.9) 



We can then consider enhancements through various fields remaining light, taking the 
mixing into account. 



5.1 Charginos 

The chargino mass matrix is expanded by new charged states from the triplet: in the 
basis iT+,W+,H+)/{T-,W-,H;l) it is 



/ Mt+^ m2D+92VT 



\tvcp 

Mch = m2D - g2VT M2 g2VSp/y/2 

With f T ~ and defining Mt = Mt + ^^^^ this has determinant 



(5.10) 



det(M, 



Ch) 



4(m^o - M2MT)fi + 2\/2g2Arm2DV^C2/3 + (A^M2 - 2g^MT) 



V^S2i3 

(5.11) 



Since the loop function A1/2 varies very little between the lower bound on chargino 
masses from LEP (105 GeV, or 92 GeV with caviats) and infinite mass, it is very well 
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approximated by 4/3 over all cases of interest. Then we can well approximate 



4 
3 



5ii5/ilogdet rrif + 52iC?h log det mj + Ssidsj^ logdetmj 



3 det Mch 



Sii ( 2V2g2XTmD2C2p + {2M2\l - giMT>2^ 



-2'S'2i f 2\/2g2\TmD2S2i3 - {2M2\t - glM^)s2i3 
--^53i(^A5(m|)2 - MsMt) + \sT{\gWs2p - M2fl) 
In the limit M2 = Mt = Xst = 0, this simplifies to 



(5.12) 



^Charginos ^ ^ \ 

■^■^ 3 y/2mD2iJ' - g2v'^>^TC2(3 _ 



2g2XTv'^{-C2i3Sii + 82(5821) + \svmD2Su 



(5.13) 



The scenario of light charginos is particularly appropriate for the MSSM in dis- 
guise; since the adjoint scalars are all massive the main phenomenological difference 
with the MSSM will be the charginos and neutralinos. Hence this represents one 
useful limit of this formula, where there is negligible mixing between the Higgs states 
(i.e. Sii ~ 1, Sii ~ V i 7^ 1). Then it is clear that for an appreciably large coupling 
At and large C2/3 the 77 channel can be enhanced without affecting the other channels. 
We performed at scan at Xt = l,tan/3 = 50 varying /i, and plot the contours of the 
Higgs to diphoton branching ratio in figure [2j showing also the effect on Ap. This 
analysis is similar in spirit to models adding a triplet to the MSSM [54j, except instead 
of including a Majorana Wino and supersymmetric triplet mass we have included a 
Dirac Wino mass. 



An alternative application of formula (5.13) is to allow appreciable Higgs mixing 



but take the large m£)2 limit, leaving light Higgsinos. We then find 



5ii -0.28 cot -0.15 — S31 



(5.14) 



Hence the diphoton production rate can be enhanced for suitably large vXs/fJ- and 
^ai. This is applicable for the NMSSM too, and is particularly interesting because 
by varying the Higgs mixing terms we can simultaneously enhance n^^, decrease the 
IJ,f£ and Hrf while maintaining fiww and nzz roughly unchanged if so desired. This 



is similar to singlet extensions of the MSSM without Dirac gauginos 55 and can be 



easily understood from equations ( 5.5 ) and ( 5.7 ) : by having a small positive admixture 



of H we enhance the coupling of the Higgs to tops, and hence to gluons; by reducing 
^ii we reduce the coupling to bottoms and taus, which reduces the total width of the 
Higgs - both of these can compensate for the reduction in coupling to Ws. In figure [3] 
we show how fihb and affected by the mixing. The above effect of the mixing 

and charginos on is then illustrated in figures |4] and [5j 
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Figure 2: Results of a scan over a subspace of the "MSSM in disguise" with At = 0.7, 
tan /3 = 50, varying mD2 and /i in order to vary the hghtest chargino mass which is restricted 
to be Dirac Wino-hke (above the contour the hghtest chargino is Higgsino-hke). The stop 
and sbottom masses were taken to be equal and varied in order to keep the Higgs mass 
at 125 ± 4 GeV (within reasonable accuracy of the results of the SARAH-produced SPheno 
code, which cannot include two- loop effects due to the Dirac gluino). Other non-zero soft 
parameters were: mm = 200 GeV, mD3 = 1000 GeV, As = 0.1; the slepton and first two 
generations of squark masses squared were 4 x 10^ (GeV)^; the singlet and triplet scalar 
masses were approximately 5000 GeV; Bfi was varied over [300,5 x 10^] (GeV)^. Left: 
contours of /x^^ (GeV). Right: /i^^ vs Ap; orange points have the lightest chargino mass 
greater than 105 GeV, blue ones smaller. 
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Figure 3: Contours of fibb or equivalently /i^r with tan/3 = 1.2. The dashed hues represent 
f^ww = 1-3 (red) , 1 (black) and 0.7 (green). 




Figure 4: Contours of /i^^ for \sv/m^+ = —2 (left) and —3 (right) with tan/3 = 1.2. 
The black dashed lines are contours of ^ww = 1.0 ± 0.3, and the red contour is a rough 
95% confidence-level preferred region according to the data given in appendix [C] (excluding 
electroweak precision data, which is dependent on the spectra of any other light particles). 
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Figure 5: Contours of fj,^^ for Xsv/m^+ = —3, tan/3 = 3 (left) and 50 (right). The 
black dashed lines are contours of ^ww = 1-0 ± 0.3, and the red contour is a rough 95% 
confidence-level preferred region according to the data given in appendix [C] (excluding 
electroweak precision data, which is dependent on the spectra of any other light particles). 



5.2 Charged Higgs 

The charged Higgs fields can also contribute to the diphoton rate in these models. 
The full charged Higgs matrix involves not just the of the MSSM, but also the 
charged triplet scalars; however, from the p-parameter constraint we know that the 
triplet scalars must be heavy (see section |4]), so we shall neglect their contribution. 
Hence we can approximate 



^Charged Higgs ^ __ _ _SuUy- + 2A| - 14A| + {g^ + gl - 2(A| + A^.)) COS 4/3 

L \ 

+ S2i{gl + gl - 2(A| + A^)) sin 4/3 

8-531 / ^^ . 1 1 . 

(5.15) 

A large contribution to the diphoton rate from charged Higgs loops can then arise 
when mixing between the lightest Higgs and the singlet is substantial. Note that light 
charged Higgs fields in the limit of very small tan (3 also often demand light stops and 
charginos to cancel large contributions to 6 — )• 57. 
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5.3 Stops and staus 



Of the squarks in the theory, the stops and staus - having the strongest couphngs to 
the hght Higgs - have been studied in the MSSM as candidates to modify the Higgs to 



diphoton rate 156 57 . In models with Dirac gauginos, the new D-term contributions 



to the potential modify the squark masses; for the stop and stau their mass matrices 
are given in appendix equations (A. 19) and (A. 20). Neglecting the normal D-term 



components, vt and the A-terms gives approximately 



Mi 



1 



m, 



1 



V2 



VyrHSi3 



m: 



+ 2gYmDivs + v'^\y1cj^ 



TR 



(5.16) 



V -^vyrHS/s 

It is straightforward to derive expressions for the the couplings of the stops and 
staus to the Higgs eigenstates, but since the full expressions are lengthy we give here 
simplified formulae neglecting subleading terms proportional to Adz and setting the 
jd-terms and vt to zero: 
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^3i{ -^vXsfi cot j3 + 



ruDiMzsw ( 
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m|^ + ml^ + ) tan f3 



(5.18) 



S: 
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--vXsfltan (3 



mpiMzsw / 

,2 



2ml 



2 I 2 

m^ + m^ 



.\/2""'" " ^ ml \ 

Here m^. , m^^ are the masses of the stop and stau eigenstates respectively. The Dirac 
mass enters into the above by shifting the mases (e.g. m?^ + m?^ = mg + + 2m^ — 
gyT^DiVs) but clearly only plays a significant role in enhancing the diphoton channel 
if there is large mixing of the lightest Higgs with the singlet scalar. 



6 MSSM without /i-term 

In the presence of Dirac gaugino masses it is possible to remove the Higgsinos' mass 
term from the superpotential of the MSSM [7]. This is another way of curing the 
intrinsic /i-problem of the MSSM. Furthermore, an approximate U{l)ii symmetry 
naturally guarantees that tan/3 is large, explaining the top/bottom quark mass hier- 
archy. In contrast to its appealing theoretical aspact, the ^SSM is under substantial 
pressure from experimental data. First, LEP put a lower limit on the mass of the 
lightest chargino of 94 GeV |58| . The chargino mass eq. ( |5.10 ) reads in this limit 



m2D Xrvcp \ 

Mch =[ m2D g2VSp/^/2 (6.1) 

-Xrvsp g2VCi3/V2 / 
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It is known that it is possible to fulfill this bound by a careful choice of At and m2D' 
a large value of At as well as m2D around 107 GeV is needed to maximize the mass of 
the lightest chargino, see Fig. [6] The highest mass which can be reached at tree-level 
is about 110 GeV. This is also not improved at the one- loop level because the loop 
corrections due to the heavy triplets even tend to decrease the mass. Hence, the 
highest mass we could find in our scans calculating the full one-loop spectrum was 
103 GeV. 




-1.0 -0.5 05 1.0 -1500 -1000 -500 500 1000 1500 

At m2D [GeV] 



Figure 6: Mass of the lightest charginos as function of At (left) and m2D (right). The 
others parameters have been chosen in the ranges: tan/3 = [30,70], Ag = [—0.2,0.2], 
X^^ = K = [-1.7,1.7], Bs = Bt = [-10^10^] GeV, 5^ = [-100,100] GeV, vt = vs = 
[— 1, 1] GeV, miT) = [—1.5, 1.5] TeV, Ttop = [—1-5, 1.5] TeV. The sfermion sector is fixed by 
^lii = = mlii = 5-106 GeV^ ml ,, = mf ., = 5 • 10^ GeV^ (^ = 1, 2), 33 = 33 = 
33 = 106 GeV^ m2 33 = ^^^33 = 10^ GeV^. 

A more severe problem is that the large values of At lead often to a huge contribu- 
tion to Ap. The reason is that this coupling breaks the custodial SU(2)l symmetry 
present in the SM and MSSM: especially, the chargino and neutralino-loops contribute 
differently to the W and Z self-energies. The large impact was already pointed out 
in ref. [T] using an approximation for the resulting contributions to the T-parameter. 
We repeat this analysis with the full numerical evaluation of Ap. We find a strong 
correlation in this model between the mass of the lightest chargino and the smallest 
possible value of Ap as depicted in Fig. [7j 

Even if the full calculations leads to somewhat smaller values of Ap than the 
approximate one for At > 1, all points which fulfill the limit of m-+ > 94 GeV 
suffer from a large Ap of at least 0.003. Many points are even above 0.01. Note that 
Ap is very quickly increasing with the chargino mass. Therefore, demanding Ap < 
0.0008 would rule out all points with chargino masses above 20 GeV. Furthermore, 
it can also be seen that in general the points with a Higgs mass between 122 and 
128 GeV lead independently from the chargino mass to Ap > 0.001. The reason 
is that a sizable contribution of At to the tree-level Higgs mass is needed. This 
could, of course, be circumvented to some extent by allowing for even larger values 
of Ttop- However, this is in contradiction to the approximate f7(l)jj-symmetry which 
suppresses in general the trilinear terms. Therefore, if we restrict ourself to moderate 
values of the squared squark mass parameters and trilinear soft-term, this model is 
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Figure 7: Left: 6p = log^Q Ap versus the lightest chargino mass. The orange points indicate 
a mass of the hghtest Higgs in the range of 122-128 GeV. Right: 6p as function of At. The 
green hne shows the resuh using the approximative formula given in ref. (?]. The input 
values are the same as for Fig. |6} 

always in conflict with Ap: even if it would be somehow possible to find kinematical 
configurations to significantly reduce the LEP limits on the chargino mass, the now 
existing bound on the Higgs mass still predicts too large values of Ap. 



7 Dynamical fi models 



Although the "MSSM without fi term" may be severely challenged, by allowing a 
substantial expectation value for the singlet the various problems can be cured. The 
Higgs potential will always lead to a non-zero value for vs as can be seen from the 



minimisation condition (A. 16), but in order for this to be significant we can allow a 
non-zero negative value for Bs and/or a non-zero tadpole term ts- Both of these are 
generically present and do not break R-symmetry. In this scenario, as in the MSSM 
without p, term, we take the only significant source of R-symmetry breaking to be a 
Bp term. 

This scenario is particularly interesting from the perspective of Higgs mixing, 
since the singlet adjoint scalar will typically be light - we see from the minimisation 



conditions that the parameter rhg in the tree-level mass-squared matrix (A. 12) is 



ml = - ^'^ + ^0 (7.1) 



and, expecting from naturalness and RGE running 34 ~ ~ v^, vs ~ v, frig ~ v'^. 
There may thus be substantial mixing between the singlet and original "/i" eigenstate; 
the size of Bp term then controls the amount of "ff" in the lightest Higgs mass 
eigenstate. 

Moreover, the singlet couples to the gauginos via the coupling Xs, which, if m£)2 
is not small, will be predominantly Higgsino-like. This then offers the possibility of 



realising the scenario considered in section 5.1 We have therefore conducted a scan 



over a portion of the parameter space of these models, concentrating on models with a 
small component of mixing between h and H but substantial and S^i copmonents. 
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using the SPheno code produced by SARAH, tan /3 was taken to be 1.5 and A5 was varied 
from a negative initial value in order to fix the Higgs mass at 125 ± 4 GeV - recall 
that this is a rather conservative error range. The other parameters varied were 
niDi e [-800,800] Ge\,vs G [130,430] GeV, S// G [312,90312] (GeV)^, while the 
non-zero fixed soft parameters were = 0.021, iJ^ = —5 x 10^ (GeV)^,t5' = —1.5 x 
10^ GeV^m2, = 2.5 x 10^ {Ge\f,m% = 9 x 10^ {G(Nf,mD2 = 600 GeV,mD3 = 
1000 GeV. The slepton and first two generations of squark masses squared were 
4x 10^ (GeV)^ while the third generation squark masses squared were 1.5x 10^ (GeV)^. 

Figure [7] shows the nww and values versus ^^-^ [.[^zz and /ibf, being, as per 
the approximate formulae, almost identical to [iww and respectively) while also 
giving the values of Ap and BR(6 — )• 57). As can be seen from the plots, there are 
many experimentally viable model points. To further elucidate the comparison with 
our predicted scenario, a plot of the mixing parameters and S'21 with //^^ revealed 
by the colour of the points is shown in figure [9} 

It is worthwhile to pick out one example point; this has A5 = — 1.3, Bfi = 
4.18 X 10^ (GeV)^'^;s = 162 GeV and ruDi = -152 GeV which leads to light neutral 
Higgs masses of 122, 227, 394 GeV, a light pseudoscalar Higgs mass of 285, a light 
charged Higgs mass of 210 GeV, neutralinos of masses 65, 156, 220, 315, 649, 651 GeV 
and charginos of masses 144, 647, 652 GeV. The mixing data is = 0.89, S'21 = 
0.08, Ssi = —0.44 which results in = 1.6,fiww = 1-14, ^zz = 0.97, /if,b = 
0.81, = 0.86 and Ap = 7 x 10-^ BR(6 ^ 57) = 3.6 x lO""^. 



Models with light stops 

In the previous scan, we held the third generation masses to be heavy to be above 
search bounds and to diminish their contribution to Ap whilst still remaining nat- 
ural. However, it is also straightforward to find models of the above class that 
have light stops, which would be natural even for the MSSM but also interesting 
for LHC searches. Taking the same fixed values as above except now with third 
generation soft masses at (500 GeV)^,m£)2 = ITeV, tan/3 = 2, one example has 
As = -0.96, 5p = 3.5 X 10^ (GeV)2,t;5 = 193 GeV and tudi = -294 GeV which 
leads to vt = 0.46 GeV, hght neutral Higgs masses of 122, 256, 360 GeV, a light 
pseudoscalar Higgs mass of 290, a light charged Higgs mass of 254 GeV, neutralinos 
of masses 115,138,319,343 GeV and charginos of masses 129,1055,1059 GeV. The 
mixing data is 5ii = 0.97, 52i = 0.03, S^i = —0.23 which results in fi^^ = 1.4, fiww = 
1.2, fizz = 1-0, Pbb = 1-0, Hrr = 1-08 and Ap = 4.8 x 10"^, BR(6 57) = 3.4 x lO"'^. 

Another example with more mixing has A5 = —1.14, Bfj, = 4.6 x 10^ (GeV)^, ^5 = 
178 GeV and moi = -283 GeV which leads to vt = 0.47 GeV, hght neutral 
Higgs masses of 121,290,380 GeV, a light pseudoscalar Higgs mass of 330, a light 
charged Higgs mass of 283 GeV, light neutralinos of masses 115,151,318,354 GeV 
and charginos of masses 140 GeV and a TeV. The mixing data is Sn = 0.95, S21 = 
0.07, S31 = —0.30 which results in fi-y^ = 1.6, pww = 1-3, pzz = 1-1, pbb = 0.95, prr = 
1.0 and Ap = 7 X 10-^,BR(6 57) = 3.4 x 10"^. 



19 




Figure 8: Production cross-section times branching ratio for a scan of a subspace of the 
"dynamical /i" scenario, with parameters chosen to enhance Higgs mixing. In the scan, tan (3 
was taken to be 1.5 and \s was varied in order to fix the Higgs mass at 125 ± 4 GeV with 
inital value being negative. The other parameters varied were mjji e [—800, 800] GeV, vs G 
[130,430] GeV,5/i G [312,90312] (GeV)^, while the non-zero fixed soft parameters were 
At = 0.021,5s = -5 x 10^ {GeVf.ts = -1.5 x 10^ GeV^ = 2.5 x 10^ (GeV)2,m^ = 
9 X 10^ (GeV)^,mD2 = 600 GeV, m^s = 1000 GeV. The slepton and first two generations 
of squark masses squared were 4 x 10^ (GeV)^ while the third generation squark masses 
squared were 1.5 x 10^ (GeV)^. The points shown in orange pass all experimental limits 
and furthermore lie within a crude 95% confidence limit via based on the current Higgs 
and electroweak precision data as described in appendix [Cj 
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Figure 9: Left: subset of results of the "dynamical /i" scan shown in figure [7] plotted in 
terms of the mixing matrix parameters Sn and 5*21 where all point pass experimental limits 
(except Higgs branching ratios). The point colours denote the /i^^ values, varying from 
dark blue for /i^^ < 1 to red with 1.33 < fi-y-y < 1.67 to green for > 2. The solid contour 
lines show fibb = 0.6, 1.0 (thick red, thick blue respectively) and Sf^ + S21 = 1 (thin blue) 
while the dashed black contours show fiww = 1-0 ± 0.3. Right: shown for comparison, this 
is a plot of the same form as figure [4] but with tan /3 = 1.5, Xv/nif = —2. 
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A final example has tan/3 = 1.5, ts = 1.2 x lO'^ {GeYf,Xs = -lAA,Bfi = 4.7 x 
10^ {GeY)^,vs = 144 GeV and niDi = -247 GeV which leads to vt = 0.3 GeV, 
light neutral Higgs masses of 125.3, 218, 389 GeV, a light pseudoscalar Higgs mass 
of 298 GeV, a light charged Higgs mass of 206 GeV, light neutralinos of masses 
89, 153, 301, 368 GeV and charginos of masses 141 GeV and a TeV. The mixing data 
is 5ii = 0.88,521 = 0.09, /Sai = —0.47 which results in fi^-y = 1.8, ^ww = 1-2, /xzz = 
1.0, fibb = 0.80, firr = 0.86 and Ap = 8 x 10''^, BR(6 57) = 3.5 x 10^1 

8 Conclusions 

Dirac gaugino models are gaining increased interest as non-minimal supersymmetric 
standard models with enhanced naturalness compared to the MSSM and an enhanced 
Higgs mass that can also relax bounds on direct superpartner searches. With the latest 
update to the SARAH package, it is now possible to study such models quantitatively 
using modern numerical tools, and this work is a first step in exploring phenomeno- 
logically the low-energy parameter space. We have discussed the properties of three 
different Dirac gaugino scenarios that are subclasses of the minimal Dirac gaugino 
extension of the (N)MSSM: the "MSSM in disguise," the "MSSM without ^ term" 
of [jl and a new scenario involving a dynamical /i term. While the first of these is 
phenomenologically very similar to the MSSM with higher-dimensional operators, we 
found that the second is unfortunately severely challenged by the current data. The 
third scenario, on the other hand, can be particularly natural and also has many 
characteristics appropriate to allow Higgs mixing and thus modifications of the Higgs 
production and decay rates; in particular, it is possible for example to enhance the 
diphoton signal, suppress the bottom and tau signals, while leaving the Z and W 
channels roughly the same as the Standard Model. We have performed a first exami- 
nation of its parameter space but clearly it would be interesting to examine it further, 
particularly as new Higgs data becomes available. 

There are now many interesting directions for future work. One will be to compare 
specific models directly with collider data, particularly in the context of models with 
light stops. In addition, it would be interesting to see how embedding the models we 
have discussed in particular high-energy completions affects the discussion of natural- 
ness. Furthermore, constraints due to dark matter (assuming a thermal history of the 
universe or otherwise) can now also be applied. On the technical side, to further refine 
the precision of the Higgs mass, the leading two-loop corrections involving the Dirac 
gluinos should now be calculated. This work is therefore one step on the increasingly 
attractive path of bringing the phenomenology of Dirac gauginos closer to the level 
of understanding of the (N)MSSM. 
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A Tree-level parameters of the model 



In this section we summarise the tree-level parameters of the model; see also 21 



Here we add the tadpole term for the singlet and expressions for the stop and stau 
mass matrices including the new Dirac gaugino D-term corrections. 



A.l Higgs potential 

It will be useful to introduce the following effective mass parameters: 
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(A.l) 



A. 1.1 Equations of motion for the CP-even neutral fields 

The scalar potential for the CP-even neutral fields is given by: 
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where the effective masses for the real parts of the S and T fields read: 
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(A.2) 



M| + ml + ^m{jj + Bs, rfiTR = + + 4m2£, + Bt (A.3) 



There is no restriction on the sign of the different mass parameters m| and Bs at 
this stage. 

The imaginary parts of the fields have been dropped as their vevs are vanishing due 
to the assumed CP conservation [13] . The coefficients of the corresponding quadratic 
terms: 



Ml + ml- Bs, 



fhj'i = + m\ — Bt 



(A.4) 



do not, in contrast to the CP-even partners, receive contributions from D-terms pro- 
portional to the Dirac masses. 

As is customarily done for the (N)MSSM, the minimization of the scalar potential 
allows here also to express jl and Bfj, as a, function of the other parameters: 



and 



P Ml 
u H 



m 



+ 



{gm2DVT - g'miDVs) 



Ml 



'^ = 212^ + ml + _ ^s+^^22 

S2I3 1^ Ha 2 P 



(A.5) 
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The new equations are 
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We can use these to solve for the masses in terms of the vevs. However, since 
the vev of T contributes to the W boson mass, the electroweak precision data give 
important bounds on the parameters of the model. For instance, using Ap = (4.2 it 
2.7) X 10^"^ I59lj63], we require: 
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< 1 X 10"'^ (95%) 



(A.9) 



which is satisfied for ^ 4 GeV. For large triplet masses we have 
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A. 1.2 Masses of the CP even neutral scalars 

Introducting the notation 
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the mass matrix for the CP even scalars in the basis {h, H, Sji, T|^} is: 

/ M^ + Ahsli^ AhS2(3C2(3 Ahs Aht \ 

AhS2i3C2i3 M\ - Ahslp Ahs Am 



\ 



Ahs 
^ht 



Ahs 



mi. 



XsXt- 



XsXt\ 



m 



T 



where we have defined: 
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which vanishes when A5 and Xt take their N = 2 values [o]. We denote non-diagonal 
elements describing the mixing of Sr and T]j states with the light Higgs h: 
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while 

^Hs = 9 miDVS2p - \s ^ C2p, 
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(A.15) 



stand for the corresponding mixing with heavier Higgs, H. 

Let us work with Mg = Mt = As = A^ = At = 0. Rewriting the vs equation we 
have 
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Then using the minimisation conditions we can write 

^hs =v[vsXs{Xs - HS213) - g'miDC20 + -/2Xsn + AsAt^t] 
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(A.17) 



Clearly A^g < M| is necessary to prevent a see-saw reduction in the lightest Higgs 
mass. We have 
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A. 2 Squark masses 

The squark masses are modified by the Dirac mass terms via the D-term contribution. 
We give here the expressions for the mass matrices for stops, sbottoms and staus: 

1 12 
{Ml)n =mQ + mf + -gymDivs + 92mD2VT + ^zi^ ~ ^^w)c2p 

{Mf)2i =mt{Atsp - ficp) 

4 2 

{M})22 =mf^ + mf - -gymDivs + -M|s^C2^ (A.19) 
The entries of the sbottom mass matrix read 
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The entries of the stau mass matrix read 



{Mf)ii =m\ + ml- gyrriDivs - g2mD2VT + M|( 
{Mf)22 =ml^ + ml + 2gYmDivs - M|s^C2^ 



Sw)c2l3 



(A.21) 



B One-loop effective potential 

The standard general form of the Higgs potential up to quartic order is 
Veff = {ml^ + i^^)\Hu? + {m\ + iJi^)\Hd\^ - [mj^Hu ■ + h.c] 
+1 [\i9' + 9") + Ai] {\H£f + 1 + g'^) + X,] 
Hlig' - g") + A3] + [-l9' + A4] {H, ■ H^){Hl . HI) 
+ {^{Hd ■ Huf + [Aeli^dP + \7\Hu\''] (Hd ■ Hu) + h.c), (B.l) 

In Dirac gaugino models where we integrate out the heavy singlet and triplet scalars, 
we find A3 and A4 have a tree level contribution, so we write A3 = 2A^ + A3,A4 = 
A| — Ay + A4 where A3, A4 are the loop corrections to the potential. Then we find at 
one loop A4 = A5 = Ae = A7 = 0, with the remaining contributions from the singlet 
and triplet scalars given by 



1287r2Ai =4A|log 



~ 2 -2 

msjv 



+ 4A|log 



~ 2 —2 

msRV 



+ 8 ^ — rh^j + nVgj 



+ 4A^log 
-1 



~ 2 -2 



+ 4A^log 



m^^ + m|^)log 



~ 2 -2 



2 —2 
UIsrV 



~ 2 —2 

rurpjv 



+ ifiTirhTRlog 



~ 2 —2 

mrpjv 



~ 2 -2 
UItrV 



+ m^jfn^Rlog 



rriTRV 



+ 4( - + fn^j^ ( ~ ( ~ "'^TR + °i^T?j ( ~ 2A|. + 



+ ^2 {zm^j — "mrjij — 4^2'^r/Af^ + ^^S ' 

+ ^2 ^ - 'irfiTR + mf./) A^ + AglrhTR\\ - gtruTR^ log 



~ 2 -2 

rrij'jV 



~2 -2 

UItrV 



) (B.2) 
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1287r^A2 =4A^log 
+ 



~ 2 -2 

mgjv 



+ 4A|log 
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-1 



~ 2 I ~ 2 



+ 4A^log 

-1 
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TR 
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\2 \2 
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+ 4 
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-1 



~ 2 I ~ 2 



2A|. + 5i 
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~ 2 -2 
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(B.3) 
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(B.4) 



Note that there are also important contributions from the stops and sbottoms, 
which at this level are identical to those in the MSSM (hence we separated them 



in equation (3.1)) except with the soft masses shifted by the Dirac gaugino D-term 



contributions as in equations (A. 19) and (A.20). 



C Experimental data 



In this appendix we give the experimental data, current at time of writing, used in 
the text for crude 95% confidence level limits. The Higgs mass reported by CMS [2] 
and ATLAS [l] is: 



mu = 125.3 ± 0.4 ± 0.5 GeV (CMS), 126.0 ± 0.4 ± 0.4 GeV (ATLAS). (C.l) 

In table [1] we give the latest reported values of na. 
The limit for Ap reported in 60 -63] is 



Ap= (4.2 ±2.7) X 10" 



(C.2) 
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CMS 


ATLAS 


Tevatron 




1.6 ±0.4 
0.64 ±0.57(7 TeV) 
0.79 ± 0.56(8 TeV) 
u.oo HZ u.ou yi lev j 
0.98 ± 0.71 (8 TeV) 
0.59 ± 1.17 (7 TeV 
0.41 ± 0.94 (8 TeV) 
0.62 ± 1.17 (7 TeV) 
-0.72 ± 0.97 (8 TeV) 


1.8 ±0.5 
1.7 ± 1.1(7 TeV) 
1.3 ± 0.8 (8 TeV) 
u.o ic u.u yi lev J 
1.9 ± 0.7 (8 TeV) 
0.46 ±2.18 (7 TeV) 

0.45 ± 1.8 (7 TeV) 


"-"■"^-2.54 

n Q9+113 
"-'■'-'^-0.32 

^■^ ' -0.68 



Table 1: Table of production cross-section times branching ratios over standard model 
values for CMS [2], ATLAS |1] and the Tevatron [64] . 



To get bounds from b — )• ,57, we define the ratio of SUSY to SM contributions 



65-70 



^ BR(6 ^ s7)susY „^ 
BR(6^.7)SM ^ ^ ^ 

Adding to the uncertainty of the SM prediction Br(5 ^s7)sM = (3.15±0.23) • 10"^ 
an intrinsic SUSY error of 0.15 as well as the error of the experimental world average 
Bt{B X^7)exp = (3.43 ± 0.22) • 10"^ (TI], leads to the following 95% CL bound 

= [0.87, 1.31] . (C.4) 
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